Stochastic resonance on the transverse displacement of swimmers in an oscillatory shear flow 
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Self-propelled microorganisms, such as unicellular algae or bacteria, swim along their director relative to 
the fluid velocity. Under a steady shear flow the director rotates in close orbit, a periodic structure that is 
preserved under an oscillatory shear flow. If the shear flow is subjected to small fluctuations produced by small 
irregularities in the microchannel or by other swimmers nearby, the director dynamics becomes stochastic. 
Numerical integration of the swimmer motion shows that there is stochastic resonance: The displacement in the 
vorticity direction is maximized for a finite noise intensity. This transverse displacement resonance is observed 
when the displacement is coarse grained over several periods, although the director is preferentially oriented 
along the flow. The resonant noise intensity is proportional to the oscillation frequency and independent of 
the shear rate. The enhanced displacement can have effects on the transverse diffusion of swimmers and the 
rheology of the suspension. 

PACS numbers: 47.63.Gd, 05.40.Ca, 47.63.mf 



Introduction. In recent decades, interest in the dynamics 
of self-propelled organisms has increased enormously. The 
interest is twofold, from continuum mechanics to describe the 
motion of single swimmers and from statistical physics to de- 
duce collective behaviors that emerge due to their mutual in- 
teractions. Self-propelled organisms belong to what is called 
active matter, in which there is a continuous energy flux from 
some reservoir to produce motion; this energy is finally dissi- 
pated via viscosity or other similar means. From a statistical 
mechanics point of view, this energy flux puts active matter in 
out-of-equilibrium conditions. 

Bacteria and unicellular algae are a particular kind of self- 
propelled organisms. Considering their micrometer scale and 
typical propulsion velocities, they are in the low-Reynolds- 
number regime, in which inertial effects are negligible in com- 
parison with viscous ones. In this category, the bacterium 
Escherichia coli (E. coli) has been intensively studied, and 
much is known about its genetics, biological processes, and 
motility |H 0] . At low Reynolds number, a swimmer can be 
modeled as a force dipole. Depending on wether the dipole 
is tensile or contractile, swimmers are classified as pushers 
or pullers, respectively. The distribution and orientation of 
the force dipoles in the fluid have rheological effects. Indeed, 
elongated swimmers placed in a shear flow orient preferen- 
tially along the extensional direction. As a result, pushers, 
like the E.coli, reduce the fluid viscosity while pullers increase 
itirSl. 

The motion of swimmers creates agitation of the fluid. 
This agitation can be visualized by placing solid particles as 
tracers. The induced tracer motion shows anomalous diffu- 
sion and in the long-time limit the induced motion is diffu- 
sive JH §]. The induced diffusion has also been observed 
close to solid surfaces ifioll . The fluid agitation affects the 
motion of nearby swimmers as well. When several swimmers 
interact, this agitation acts as a self-induced noise f7llrl - 131. 
In this Brief Report we study the effect of flow noise (either 
self-induced noise or that from other sources) on the swim- 
mer motion under an oscillatory shear. It will be shown that 



the noise enters multiplicatively in the swimmer equation of 
motion which is itself nonlinear. 

An interesting phenomenon called stochastic resonance 
(SR), in which some response function is maximized for a 
finite noise intensity, can appear when a system is forced pe- 
riodically. It first appeared that bistability, periodic forcing 
and random forces were necessary for the onset of SR 114 1. 
However, it later became clear that SR may appear in a large 
variety of systems, including linear systems subjected to mul- 
tiplicative noise rather than to additive noise Il5lll6ll . In this 
Brief Repport we present a SR in swimmers: In an oscilla- 
tory shear flow, the displacement in the vorticity direction is 
maximized for a given noise intensity. 

Deterministic swimmer dynamics. Consider a self- 
propelled swimmer moving at low Reynolds number in a 
fluid. The swimmer propels with respect to the fluid at a 
velocity Vo pointed by the director n. If the swimmer's body 
is much smaller than the typical distance in which the fluid 
velocity changes, the Faxen correction can be neglected and 
the total velocity of the swimmer is V = Vb« + v, where v is 
the fluid velocity at the center of the swimmer 

OIH. The 

velocity gradients induce rotation of the swimmer, described 
by Jeffery's equation for the director vector lfl9ll 

h = (l-hn)[pE s + E AS }h. (1) 

Here E s and E AS are the symmetric and antisymmetric parts 
of the velocity gradient tensor E = Vv and j3 depends on the 
geometry of the swimmer. Limiting cases are /3 = 1 for a 
rod-like body and J3 = — 1 for a disk-like body. Here we are 
interested in the effects of imposing an oscillatory flow on the 
swimmer motion. 

Under experimental conditions, it is possible to place bac- 
teria in a Hele-Shaw geometry and impose an oscillatory flow. 
The resulting Poiseuille flow has a non uniform shear rate 
and, as a consequence of Eq. (Q]), the orientation dynamics of 
the bacteria depends on the vertical position in the cell. It is 
known that pushers tend to approach solid surfaces and swim 
close to them, remaining trapped by the surface for long pe- 
riods USUI. This fact allows us to simplify the analysis. 
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Close to the surface the shear rate is roughly constant and the 
same in the top and bottom surfaces. Therefore, a large ma- 
jority of the bacteria swim in the flow with the same shear 
rate. To separate the effects of walls from the effect of the 
oscillatory shear flow, we will consider a simple shear flow 
throughout the fluid. Although the equations of motion do not 
distinguish between pusher or puller swimmers, the previous 
discussion suggests that the simple shear approximation is of 
more relevance for pushers. Finally, we neglect the circular 
motion, with radius of some tens of microns, that flagellated 
swimmers develop near surfaces 112 lll23l 12411 . 

The swimmer is placed in an imposed simple oscillatory 
shear flow v°° = Y° cos(of)yi (Fig. [T]). Directions are such 
that x is the flow direction, y is the gradient direction, and z is 
the vorticity direction. 




FIG. 1. Model setup. The swimmer moves in a fluid with an im- 
posed oscillatory shear flow v°° = y°° cos (cot)yx. The dashed curve 
depicts a possible trajectory and the coarse grained displacement A 
is evaluated for lapse times of several periods. 

In the case of a stationary flow (o = 0), the temporal evo- 
lution of the director depends strongly on the value of j3 . Fig- 
ure|2]shows the evolution of the director in the cases j3 = 0.6 
and j3 = 1 . In the first case, the orbits are closed with a pe- 
riod proportional to 1 / f°, while in the second case the orbit is 
open and the director orients asymptotically to the ±x direc- 
tion. On the other hand, when the flow is oscillatory the orbits 
are always closed (periodic) and behave qualitatively similarly 
for different values of j3 as shown in Fig. [2] Based on these 
results we will analyze the case j3 = 1 in what follows. Other 
cases were studied showing qualitatively similar results lf26Tl . 

Noisy oscillatory flow. When a microscopic swimmer is 
placed in a flow, it is subjected to fluctuations of different ori- 
gin. First, there is the thermal (Brownian) force that produces 
random reorientations and thermal components on the veloc- 
ity. The effects of these have been largely studied and result in 
an effective diffusive motion J2l and smoothing of the Jeffery 
orbits fiH. 

The fluctuations on the velocity field (and velocity gradi- 
ent) in which the swimmer moves are another source of noise. 
Velocity fluctuations do not produce large effects except on an 
added diffusive motion of the swimmer. However, as it will be 
shown, fluctuations in the velocity gradient lead to a preferen- 
tial displacement of the swimmer in the vorticity direction, an 
effect that is magnified for a given noise intensity. 

In experiments, the noise in the velocity gradient can have 
several sources, three of which we mention here. First, the 
imposed oscillatory flow can deviate from a perfect sinusoidal 




FIG. 2. (Color online). Time evolution of the components of the 
swimmer director: n x (gray lines, red online), n v (light gray lines, 
green online) and n z (dark gray lines, blue online). The swimmer 
geometric parameters are j6 = 1 (top) and /3 = 0.6 (bottom). The 
main figures present the case of an imposed oscillatory shear flow 
with ft>/y°° = 2k, and the insets present the case of a steady shear 
flow with co = 0. The curves depend on the initial conditions but are 
qualitatively similar for other initial conditions. 



in an uncontrolled way. Second, the micro-channel can have 
(sub) micrometric roughness inducing velocity fluctuations in 
the Lagrangian frame of the fluid. Finally, other swimmers 
in the vicinity of the studied object create currents that are 
superimposed to the oscillatory flow. The intensities of these 
fluctuations depend on experimental conditions and swimmer 
concentration 111 ll4l3ll . 

In the presence of noise, the shear rate tensor that ap- 
pears in Jeffery's equation (HJ is E(f) = E°°(r) +E noise (f), 
where E°°(f) = Y° cos((Ot)xy corresponds to the imposed 
oscillatory shear flow and E nolse (f) takes into account the 
velocity fluctuations. It is modeled as a tensor of white 
noise components of intensity F, with the trace subtracted to 
model an incompressible flow. That is, an intermediate ten- 
sor is built with components Fji satisfying (fy (£)/*£/(/')) = 
T8jkO~ji8{t — t'). Then, the tensor E nolse has components 
Fit — Fjjdjk/3, and the resulting correlations are 



£noise _ 

(% ise (04f se (O} = r (SikSjt - StjSa/3) S(t - 1'). In sum- 
mary, the position r and director h evolve according to the 
following equations 
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The noise intensity F has units of inverse of time. It should 
be compared either with 7°° or with co to quantify if the noise 
is large or small. Considering the shear rate 7°°, the oscilla- 
tion frequency co and the noise intensity F, two dimensionless 
parameters can be varied. We chose to fix 7°° and vary F and 
CO. Four frequencies are used, ©1/7°° = ft/15, 002/7°° = ft/10, 
0*3/7°° = ft/5, and 0)4/7°° = 2ft/3, while F/co is varied in a 
wide range. We recall that the explored values of co/y are ex- 
perimentally feasible. For example, in a microfluidic device 
of cross section L z x L y = 300 x 50/im 2 , with an imposed flux 
2 = 5 nl/s, the studied frequencies scan the range co <~ 0.6 — 13 
Hz. 

In the equation for the swimmer director the noise is mul- 
tiplicative. This can lead to complex phenomena in contrast 
to the effects produced by the additive noise that represents 
the thermal noise 1 1-411 - The purpose of the present study is to 
describe the effect of varying the noise intensity. 

Equqations (f2]) and (0) are interpreted according to the 
Stratonovich calculus, and they are numerically integrated 
using Heun's predictor-corrector method 1271 12811 . Figure [3] 
shows the evolution of the components of the director vec- 
tor for different noise intensities. For small noise intensities, 
the Jeffery orbits are slightly perturbed. For long times (not 
shown in the figure) the director oscillates around the same 
direction fixed by the initial conditions. For large noise inten- 
sities the periodic structure is completely lost and the director 
performs a random motion. In the case of an intermediate 
noise intensity, the orbits preserve some periodicity, but it is 
lost after some periods. In the long-time term the swimmer 
oscillates for long periods around some orientations until it 
switches to a new orientation. The transitions take place at 
random times and the new orientations are also random. 

To quantify the degree of orientation, averages of the 
quadratic components of the director vector are considered. 
The swimmer preferentially orients along the x axis except 
for large noise intensities {F / 7°° ~ 1) when the swimmer ori- 
ents isotropically (see Supplemental Material J3]). Note that 
in the case of passive elongated fibers, it has been found that 
they preferentially orient along the z axis l29tl . However, in 
that case, the alignment is produced by excluded volume ef- 
fects that are absent in our case because we consider an iso- 
lated swimmer. 

Although the swimmer orients principally on the x axis 
(positive and negative directions), the mean displacement can 
be in a different direction as the oscillations can lead to can- 
cellations in the x direction. To subtract the effect of the rapid 
oscillations, the displacement vector A is computed for lapse 
times of 30 periods, and later divided by 30. This number of 
periods is sufficiently large to obtain a coarse-grained descrip- 
tion of the displacements, averaging over the back-and-forth 
motion induced by the Jeffery orbits and the shear flow os- 
cillations, but small enough to capture the coherent motion 
shown in Fig. [3jc). From the displacement vector, averages 
(A 2 ) and (A 2 ) are computed. The average (A 2 ) is not well de- 
fined because it depends on the streamline in which the swim- 
mer is located. 




FIG. 3. (Color online) Time evolution of the components of the 
swimmer director under a noisy oscillatory shear flow: n x (gray lines, 
red online), n y (light gray lines, green online), and n z (dark gray lines, 
blue online). The noise intensities are (a) r/y°° = 0.001, (b) T/f° = 
0.01 (resonant noise intensity), (c) same value of T for a longer time, 
and (d) T/f = 0.316. The swimmer geometric parameter is j3 = 1 
and the frequency of the oscillatory shear flow is oh. The curves 
depend on the initial conditions but are qualitatively similar for other 
initial conditions. 



The average squared displacements are shown in Fig. [4] as 
a function of the noise intensity. The displacement in the vor- 
ticity direction-(A?)-shows a maximum for a small but finite 
noise intensity. The swimmers show an enhancement of the 
transverse motion even though the director is mainly oriented 
along the x direction. Figure H shows the resonant displace- 
ment normalized with the distance traveled in one period, pre- 
senting a weak increase with the oscillation frequency. The 
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squared displacements in the y direction are smaller than those 
in the z direction and show no maximum. 

The observed phenomenon is a stochastic resonance in 
which the response (transverse displacement) is maximized 
for a finite noise intensity 

BEES. 

Smaller noises lead 

to slightly perturbed Jefferry's orbits and large noises produce 
a complete isotropic response. At the resonant noise inten- 
sity the trajectory has the appearance shown in Fig. 0c). It 
should be mentioned that the resonant noise intensity is small, 
p es <c 7°°, 0). This is responsible for the large transition times 
observed in Fig. [2c). The resonance curves for different fre- 
quencies collapse when plotted against T/co and rescaled to 
their maximum value. The resonance is wide and rather flat, 
making it difficult to identify the resonance noise intensity r les 
with precision. It lays in the range r res = (8 — 30) x 10~ 2 a>. 

Stochastic resonance was observed in linear systems (i.e., 
without bistability), subjected to multiplicative noise as long 
as the noise had some finite correlation time IU5U16I1 . In the 



present case, the noise correlation time is zero but two aspects 
could have allowed to overcome this limitation. First, the 
equations are non-linear and no simple analysis excludes SR 
in this case. Secondly, this is a coupled system of equations 
and it is known that in this case finite correlation times can de- 
velop as in the Langevin modeling of the Ornstein-Uhlenbeck 
noise l3lll . 
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FIG. 4. Average squared coarse grained displacements as a function 
of the noise intensity T. The linear vertical scale has been rescaled 
for each frequency to the maximum value of (Ar). The oscillation 
frequencies are W\ (triangles), ah (circles), ©3 (squares), and 0)4 (di- 
amonds). (A?) in open symbols and (A^) in solid symbols (only one 
frequency is shown, others are similar). Inset: Maximum squared 
coarse-grained displacement in the z direction as a function of the 
oscillation frequency, normalized with the distance traveled in one 
period. 

Perspectives. The dynamics of a single swimmer moving 
at low Reynolds number is studied in the presence of an os- 
cillatory shear flow. In microfluidic devices, it is known that 
bacteria approach and swim close to the solid surfaces. There- 
fore, the flow acting on a swimmer can be modeled as a simple 
shear. When flow fluctuations are taken into account, the re- 
sulting equation for the director vector is non-linear and the 
noise is multiplicative. When analyzing the coarse-grained 



displacement, stochastic resonance is observed. The displace- 
ment in the vorticity direction, transverse to the flow, is max- 
imized for a finite noise intensity. The resonance is wide, 
therefore it is difficult to obtain precisely the resonant noise 
intensity P es . Based on the results for four different frequen- 
cies it is found that r res °< co, with a small proportionality con- 
stant. 

The observed stochastic resonance implies that if the trans- 
verse motion were diffusive, the transverse diffusivity and 
mixing could be maximized by varying the noise intensity. 
Also, the preferential orientation of the director along the flow 
can have rheological effects on the frequency-dependent vis- 
cosities, an analysis that is being performed and will be pub- 
lished elsewhere. 

In the model, we have not considered the rapid reorientation 
of bacteria (tumbling). The effect of tumbling would be to 
render the motion more isotropic, thus reducing the resonance 
amplitude. In addition, the hydrodynamic interactions with 
surfaces could affect the resonance. 

Finally, a more detailed analysis of the stochastic process is 
necessary to identify the key ingredients that produce the res- 
onance. It is, however, complex, as there are several degrees 
of freedom and two intrinsic frequencies, (0 and y. A simpler 
toy model could provide insight. 
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